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Abstract. 

We consider properties of the operators D{r^ M) = a^{a^a)^ (which we call generalized 
Laguerre-type derivatives), with r = 1,2,... , M = 0,1,... , where a and are 
boson annihilation and creation operators respectively, satisfying [a^a^] = 1. We 
obtain explicit formulas for the normally ordered form of arbitrary Taylor-expandable 
functions of D{r^ M) with the help of an operator relation which generalizes the 
Dobinski formula. Coherent state expectation values of certain operator functions 
of D{r^ M) turn out to be generating functions of combinatorial numbers. In many 
cases the corresponding combinatorial structures can be explicitly identified. 
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1. Introduction 

Among many ways of generalizing the ordinary derivative the notion of the so-caUed 
Laguerre derivative [1] seems to be particularly fruitful. The idea is to extend the 
operator ^ to a simple homogeneous counterpart which we define as in [2], [3] (note 
that here we omit the factor (—1) present in these references): 

= ±x-. (1) 

^ dx dx 

In Refs.[l],[2],[3] many important consequences of the replacement ^ ^ i^^; in the 
integral transform methods and in the operational calculus were investigated. The 
link between and the Laguerre polynomials becomes clear if one notices the 
operational relation (see Eq.(5) of Ref.[3]) which is easy to get using the amusing identity 

d: = {±T--{±T [4] 

e^^-x^ = n!2/^L.(--), (2) 

y 

where L^^z) are Laguerre polynomials. This justifies a posteriori the name Laguerre 
derivative for D^- Using Eq. (2) we may obtain the action of e^^^ on various functions, 
using the diflFerent generating functions of Laguerre polynomials listed in Section 5.11 of 
[5]. In particular, using the well known ordinary generating function of I/^(x) (formula 
5.11.2.1 for a = of [5] ) one obtains [6] 

valid for |6A| < 1 [6]. Analogously, using the formula 5.11.2.6 of [5], one gets 

e^^^Fi {[hi [llx) = ,F^ ([5], [1], ^) , (4) 

with iFi the hypergeometric function:!: which for many values of b specializes to 
elementary or known special functions. Note that for both these examples the action 
of e^^^ results in a substitution and a prefactor which is reminiscent of the so-called 
Sheffer-type operators [8]. 

We now employ the operational equivalence 
d 

[— ,x] = 1 < — > [a,a^] = 1, (5) 

where a, a'^ are boson annihilation and creation operators respectively and rewrite 
Dx < — > D as 

D = aa^a. (6) 

By going one step further we extend Eq.(6) by defining the generalized Laguerre 
derivative D{r^M) as 

D(r,M) = a'-(ata)^^D,(r,M) = (^r(a;^r , r = l,2,... , M = 0, 1, ... . (7) 

:|: We use a convenient and self-explanatory notation for the hypergeometric functions of type pFq\ 
pFg([List of p upper parameters] , [List of q lower parameters], x). 
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These operators are the object of our present study. Although the equivalence in Eq.(7) 
between D{r^M) and Dx{r^M) is formal since the domains of a, a'^ and x are 
different, we shall show that it provides one with an effective calculational tool. 

Since (a'^a)^ conserves the number of bosons, the operators D{r^M) act as 
monomials in boson operators which annihilate r bosons. Recent experiments in 
quantum optics have shown how one may produce quantum states with specified 
numbers of photons. This in turn raises the interesting possibility of producing exotic 
coherent states; that is, states other than the standard ones which satisfy a\z) = z\z) 
[9]. The current work introduces operators whose eigenstates may be used to model new 
coherent states which have many of the features of the standard ones, and still permit 
explicit analytic description. The explicit forms of these new generalized coherent states 
can be used to evaluate relevant physical parameters, such as the photon distribution 
and the Mandel parameter, squeezing factors and signal-to-noise ratio, etc. 

Much theoretical work has been devoted to the description of nonstandard coherent 
states; for example, the so-called nonlinear coherent states [10], multiphoton coherent 
states [11] and g-deformed coherent states [12]. The structure embodied in definition 
Eq.(7) is a special case of the extension of boson operators proposed in the construction 
of nonlinear coherent states [10]. In this latter reference one defines the generalized 
boson annihilator h by 

h = af{a^a) (8) 

choosing the /(x) that most suits the problem in question. Evidently for this 
identification r = 1 and /(x) = . In this case the commutator is equal to 

[D{1, M), D\l, M)] = {a^a + 1)^^+^ - {a^af^^\ (9) 

This emphasizes the fact that although D{1^M) and D^{1^M) annihilate and create 
one boson, respectively, they are not canonical boson operators (unless M = 0). Eq.(9) 
is a special case of 

/r+l 

[D{r, M), D\r, M)] = {a^a + rf^ ( \a{r + 1, k)\{a^a)^-^ 




(a^a^ > Ja(r,fc)|(ata)M, (10) 



where the a(r, k) are Stirling numbers of the first kind [13]. 

Eq.(lO) was obtained by using the following two equations 

a'{a^Y = n (a^a + p) = 
p=i 

= '^^\(T{r + l,k)\{a^a)^-\ 
k=i 

The first part of Eq.(ll) is readily proved by induction. To prove the second part of 
Eq.(ll) we use the generating function for |cr(r + 1, in the form [14] 

EW{r + l,k)\x'+^-^ = fl{l+px) (12) 

k=l p=l 



(11) 
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from which, by substituting x = 1/n and using the first part of Eq.(ll), the second part 
of Eq.(ll) follows. 

The basic objective of this work is the investigation of arbitrary powers of D(r, M) 
which in turn will allow one to evaluate Taylor-expandable functions of D{r^M). We 
achieve our goal following recently developed methods of construction of normally 
ordered products [15], [16], [17], [18]. As we shall show, results derived in this way 
have a combinatorial flavour and lend themselves to a combinatorial interpretation. 

The paper is organized as follows. In Section 2 we introduce generalizations of the 
Stirling and Bell numbers which are well known from classical combinatorics and relate 
them to the normally ordered powers of operators D(r, M). These numbers, as shown 
in Section 3, may be explicitly found using generalized Dobihski relations. In Section 
4 we compare calculations of purely analytical origin with those based on methods of 
graph theory and give a combinatorial interpretation of our results. Examples of various 
applications of our approach are presented in Section 5 while Section 6 summarizes the 
paper. 



2. Normal ordering: Generalized Stirling and Bell Numbers 

The normally ordered form of F(a, a'''), denoted by Fj^{a^ a'^) [19] is obtained by moving 
all annihilators to the right using the canonical commutation relation of Eq.(5). It 
satisfles Fj^{a^a'^) = F{a^a'^). On the other hand the double dot operation : G{a^a'^) : 
means that we are applying the same ordering procedure but without taking account of 
the commutation relation. Conventionally the solution to the normal ordering problem 
is obtained if a function G(a, a'^) is found satisfying 

F^{a,a^) = :G{a,a^) : . (13) 

A large body of research has been recently devoted to flnding the solution of Eq.(13) 
[20]. A general approach which facilitates a combinatorial interpretation of quantum 
mechanical quantities is to use the coherent state representation. Standard coherent 
states 



oo 



.) = e-NV^5:4=|n) (14) 



n=0 V ^! 



with the number states \n) satisfying a'^a\n) = {n\n^) = 6n^n' and z complex, are 

eigenstates of the annihilation operator, i.e. a\z) = z\z). The latter eigenstate property 
shows that having solved the normal ordering problem Eq.(13) for an operator F(a, a^) 
we immediately flnd 

{z\FMia,a^)\z) = G{z,z*) . (15) 

An early observation on how to extract combinatorial content from normally ordered 
forms [21] was based on the formula e^^^^ = : e^^^^^^~^^ : [22]. It led to the identiflcation 

{z\{a^aY\z) B{n), (16) 
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where the B{n) are conventional BeU numbers described in [13]. Eq.(16) may be taken 
as a definition of the BeU numbers. For Stirhng numbers of the second kind we have 
[13], 

n 

[a^aY = ^S{n,k){a^fa^ (17) 

(which may also be used as a practical definition) in terms of which one defines the Bell 
polynomials by 

n 

B{n,x) ="^S{n,k)x\ (18) 

k=0 

We have extended and developed the coherent state methodology for operators 
other than a'^a in [15], [16] and [17]. 

After the seminal observation by Katriel [21], combinatorial methods found 
widespread application in this context [15], [16], [17], [18], [23]. We apply these methods 
to F(a, at) = [D{r, M)]^, n = 1, 2, ... . 

Formally, we write [D{r^ M)]^ in normally ordered form as 

~ Mn 

[i^(r,M)]"= ^Sl^\n,k){a^Ya^ 

.k=0 

Clearly, from Eq.(19) the integers are generalizations of the conventional 

Stirling numbers of the second kind (which are recovered for r = 0, M = 1). Analogously 
to Eq.(18) the numbers Sr^\n^ k) serve to define the generalized Bell polynomials 

Mn 

Bl^\n,x) =Y,Sl^\ri,k)x\ (20) 

k=0 

Finding the explicit form of these generalized Stirling numbers will give the normally 
ordered form of [D{r^ M)]^. We proceed to do this in the next section by use of a 
generalization of the famous Dobinski formula. 



a"^. (19) 



3. Generalized Dobinski formula 



We first write Eq.(19) in derivative form as 



Mn 



.k=0 



d 

dx 



d 

dx 



(21) 



Acting with the r.h.s. of Eq.(21) on one obtains Br^\n^ x)e^. The action of the 
l.h.s. of Eq.(21) on is obtained by acting with generalized Laguerre derivatives on 
monomials x^ 



r M p—r 



(22) 
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where p- = p{p — l)...{p — r + 1) is the faUing factorial, then extending it to the n-th 
power 

~n-l 

[D,(r,M)]"x^' = 



X 



n—rp 



YliP - rj)-{p - rj) 
ij=o 

and next summing up contributions for x^/p\ 

Tn-l 

YliP - rj)-{p - rj) 



(23) 



p=0 p=rn 



M 



J=0 



.p—rn 



p\ 



(24) 



Upon simplifying Eq.(24) leads to the Dobinski-type representation of generalized Bell 
polynomials [18], [24], [25]: 

X 

J. ' 



1=0 



(25) 



verified by direct calculation of (;2|[a''(a^a)^]"|;2). The classic Dobinski formula [13] 
corresponds to r = 0, M = 1: 

°° in I 

B{n,x) = e-J2 



1=0 



n 



(26) 



From Eq.(25) the generahzed Stirhng numbers are obtained by standard Cauchy 
multiphcation of series 



1 



j=0 



k 
j 



(-1) 



k-j 



-i M 



ir) 



i=l 



k = OA 



, J., . . . , 



Mn. (27) 



We point out that Eqs.(25) and (27) are the central results we need for further 
calculations. For practical applications it is useful to note that the generalized Stirling 
and Bell numbers, as well as generalized Bell polynomials, can be expressed through 
generalized hypergeometric functions pFq. 

Below we quote some examples of such relations. 

;-l)^(n!)^ 



M times 



M times 



• M+iFM{[-k, n + 1, n + Ij , [1, 1], 1) (28) 



B[''\n,x) = e--(n!)^ ■ mFmHu + 1, .. .,n + Ij , [l^],x) (29) 

M times M times 



The numbers B^\n) = B[^\n^l) can be shown to be related to the numbers 
Bp^p{n) (introduced in Refs.[15] and [16]) characterizing the normal order of [{a'^Ya^]^ 
by the formula 

5i^)(n)=5„,„(M + l) , (30) 
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as seen by comparing Eq.(2.6) in Ref.[16] with Eq.(25) of the present work. 



Bf'\n,x) = ^2^"e-- (n!)^ • MFM+i( [n + 1, n + Ij , 1/2], xV4)7r^/2- 

\ M times M times 

+2^ [T{n + 3/2)f ■ X ■ MFM+i([ n + 3/2, n + 3/2 ], [ 3/2, 3/2 ], a:V4) j , 

M times M+1 times / 

\n^x) = 2M+i[7^r(2/3)]^ ^ ^ ^ 

2M+i3Mn[^^,r(2/3)]^ ■ MFM+2( [n + l,. . .,n + lJ , [1^^.., 1, 1/3, 2/3], xV27)+ 



(31) 



M times M times 



+2 • 3M(n+m (r2(2/3)r(n + 4/3))'^ X ^32) 
X a; ■ MFM+2{[n + 4/3, n + 4/3], [4/3, 4/3, 2/3], a;V27) + 



M times M+1 times 

1^ ^2 



^3M(n+i) [^r(^ ^ 5/3^]M ^ ^2 . ^FM+2([ p + 5/3,. . .,n + 5/3 ], [ 5/3, . . .,5/3 , 4/3], xV27) 

M times M+1 times 

We conjecture that in general bI^^ (n, x) is a combination of r hypergeometric 
functions of type mFm+t-i of argument jr^ . 

Examples of numbers resulting from Eqs.(29)-(32) for n = 0, . . . , 6, . . . are 

M = 1 ^ (n) = 1, 2, 7, 34, 209, 1546, 13227 . . . , 
M = 2 5f ^ {n) = \, 5, 87, 2971, 163121, 12962661 . . . , 
M = 3 5f ^ (?^) = 1, 15, 1657, 513559, 326922081, 363303011071 . . . , 
M = 2 Bf (n) = 1, 10, 339, 23395, 2682076, A^l\Ylhl\ 107943795145 . . . , 
M = 3 Bf (n) = 1, 37, 9415, 7063615, 11360980081, 33040809105661, 

156151310977544887... , 
M = 3 Bf (n) = 1, 77, 39839, 62310039, 214107236041, 1358185668416501, 

14247249149298651007. . . , 
M = 4 Bf (n) = 1, 372, 1905633, 43249617004, 2805942285116705, 



(33) 



411223445534704016116, 117428972441699060660584977 . . . , 
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which are positive integers and as such admit combinatorial interpretation. The first 
two sequences in Eq.(33) may be identified as A002720 (which enumerates matching 
numbers of a perfect graph K{n^n)) and A069948, respectively, in Ref. [26]. 

We note in passing that the numbers Br^\n) are solutions of the Stieltjes moment 
problem, i.e. they are the n-th moments of positive weight functions on the positive 
half axis. This can be deduced from their Dobihski-type relations Eq.(25), whose form 
allows one to obtain the weight functions for any r and M. For the first two sequences 
in Eq.(33) the Stieltjes weights are given in [26] under their entries. 

As a second illustration of our approach we shall apply it to i^(r, 1). Note that 

[a"(ata)^]" = : Bl''\n, a^a) a'^ : , (34) 
which upon using the Dobihski relation Eq.(25) for M = 1 leads to 

The operator D{r^ 1) is of SheflFer-type viewed through hermitean conjugation (see Refs. 
[8], [27]) and Eq.(35) can also be obtained through the methods developed in Ref. [8] (see 
Appendix) . Consequently, 

'S'^exp(^^-l). (36) 
- EBi'HnA. (37) 

n=0 



where 



{z\ {a'a^ay \z) ^^(n) = \a{n + 1, p)|r"-^+^5(p - 1) . (38) 

p=i 

In Eq.(38) a{n^k) are Stirling numbers of the first kind, B{n) are conventional Bell 
numbers, and in obtaining Eq.(38) we have again used Eq.(ll). 

Using Eqs.(19) and (29) we obtain for r = 1 the following formula in a compact 
notation 

[Z)(l,M)r= ■ g-atg . ^p^^ ^ ^ 1^ ^^^y 1], gt a)a" : .(39) 

M times M times 

The last formula can be used to normally order H{\D{1^ M)) for any Taylor-expandable 
H{x). 



4. Combinatorics of normally ordered Laguerre derivatives 

In previous Sections we considered the normal ordering of Laguerre derivatives for which 
the results heavily exploited combinatorial identities stemming from the underlying 
iterative character of the problem. Indeed, the reordering of the operators a and a"^ is a 
purely combinatorial task which can be interpreted in terms of graphs [18], [28], [29] and 
analyzed by the use of combinatorial constructors [30] . Briefiy, to each operator in the 
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normally ordered form H = ^^gOf^,^ (a^)^a^ one associates a set of one- vertex graphs 
such that each vertex • carries weight ar^s and has r outgoing and s incoming lines 
whose free ends are marked with white o and gray o spots respectively. Multi-vertex 
graphs are built in a step-by-step manner by adding one vertex at each consecutive 
step and joining some of its incoming lines with some the free outgoing lines of the 
graph constructed in the previous step. Additionally, one keeps track of the history by 
labeling each vertex by the number of steps in which it was introduced. As a result, one 
obtains a set of increasingly labeled multi- vertex graphs with some free incoming and 
outgoing lines. It can be shown that the normal ordering of powers of the operator H 
can be obtained by enumeration of such structures. Namely, the coefficient of [a'^Ya^ 
in the normally ordered form of the operator is obtained by counting all possible 
graphs with n vertices • and having k white o and / gray o spots respectively. For 
illustration, we give two examples of Laguerre derivatives D{1^ 1) = aa'^a = a'^a^ -\-a and 
D(2, 1) = a^a'^a = a'^a^ + 2a^ and their graph representation leading to the solution of 
the normal ordering problem by simple enumeration (see Fig. 1). One should compare 
these "graphical results" with the explicit formulas of Eqs.(27) and (25) or the expansion 
coefficients of the generating function in Eq.(36) for r = 1,2 and M = 1. Thus, using 
Eq.(33), the coefficients multiplying the operators in Fig. (la) are the first two terms in 
B[^\n) = 2, 7, 34, 209 forn = 1, 2, . . . (A002720). Similar coefficients in Fig.(lb) are the 
first two terms in B^^\n) = 3, 16, 121, 1179 forn = 1, 2, . . . (A121629). 



5. Examples 

1. For r = M = 1, i.e. for D{1^ 1) = aa'^a one obtains (see [31] for a similar calculation): 

[D(l,ir = n!:L„(-ata):a", (40) 

where Ln{y) are Laguerre polynomials and Eq.(40) is derived from Eq.(39) and using 
the definition of Z/„(y) via the function iFi. Then 



^ ^oo 1 / Xa^a^ \ ^^^^ 

where in Eq.(41) we have used the ordinary generating function for the Laguerre 
polynomials [5]. 

Using other generating functions listed on p. 704 of Ref.[5] one can derive further 
formulas of type Eq.(41). ( In a) and b) below: A 7^ 0, p = 1, 2, ...). 

a) Formula 5.11.2.6 of [5] for a = provides the normal ordering of 

.F. (16), [II, AD(1, D) = -.j^^.F, HI. ^) : . (42) 
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which for b integer and half-integer can be written down in terms of known functions. 
Examples are: 

1 



iFi([3],[l],AD(l,l)) 



(1 - Aa)3 



exp 



1 - Aa 



1 — Aa 



iFi([|],[l],AD(l,l)) = :exp 



X 



2(1 - Aa) 



X 



(43) 



^2(1 - Aa) 

where /o(y) and l\{y) are modified Bessel functions. 

b) Similarly, we consider the formula 5.11.2.8 of [5]: 

1 



E 



n=0 



n + p 
n 



_ tx 



X 



l-t 



(1 - t)P+i 

Using Eq.(40) we obtain the normally ordered form of [AD(1, 1)]^ exp(AD(l, l))/p! : 
1 [AD(1,1)]"+P 



(44) 



{n + p)\ / t ^M ^ 



n=0 



plnl 
1 



-e Lr, 



(45) 



■) : {Xaf . 



■ (1- Aa)P+i^ 1- Aa' 

2. The normal order of the modified Bessel function of the first kind 
/o (2(AD(1, 1)Y/^) may be derived: 

L„(-ata)(Aa)" 

(46) 



/o(2(AD(l,l))V2) ^ . 



n=0 



= : e^Vo(2yA(-aata)) : = : 6^*^/0 ( 2 ^Aot^) : , 

where in the last line we have used the exponential generating function of Laguerre 
polynomials [5]. The analogous formula for Jq (2(XD(1, 1))^^^) reads 

L„(-ata)(-Aa)" 

~ ' " (47) 



Jo(2(AD(l,l))V2) ^ 



n=0 

= : e-^"/o(2VAata2) : . 

3. We quote here the eigenfunctions of Drc{r,M) with eigenvalue 1 satisfying 
Dx{r, M)E(r, M; x) = E(r, M, x), with the following r boundary conditions: 



E{r,M;0) 



dxP 



E{r,M;x) 



= 0, p=l,...,r-l 



(48) 



x=0 



which are 



E{r,M;x) = oFM+r-i([ ], [ }/r, 2/r, (r - l)/r, 1, Ijx^y^^) . (49) 

r— 1 times M times 

Useful normal ordering formulas can be obtained by applying the Dobinski relations to 
the eigenfunctions of Dx{l^M) with the argument taking operator values^ see Eq.(49), 
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i.e. E{l^M]Dx{l^M)). We briefly show the calculation, in boson notation, for 
^(l,2;AD(l,2)) =oF2([ ], [1, 1], AD(1, 2)), see Eq.(39) : 

oF,{[ ], [1, 1], AD(1, 2)) = E 7-ni[«(«^«)T = 

n=0 \^-) 

.:e--g<4;^E^^^(A<.)"- (50) 

= : e-ata g ^^^([1 + 1,1 + 1], [1, 1], Aa) : , 

and similarly 

E{l,M-\D{l,M))= oFm([ ],[M^],AD(1,M)) = 

M times 

= : e-"^" ^ ^ Mi^M([ 1 + /, 1 + ^ 1 + / ], [1, 1, ■-, IJ, Aa) : , 

M times M times 

which indicates a pattern appearing in the course of this procedure. 

Indeed, by evaluating the coherent state expectation value of Eq.(51) between 
(z = 1| . . . |z = 1) in the spirit of Eq.(16) we obtain the hypergeometric generating 
functions of the numbers (n) as then 



oo 

E 7T MFM{ [l + l,l + l,...,l + lj , [M,_J, A) 



_ 11 

'-—^ M times M times (^52) 

The hypergeometric generating function of Br^\n^x) for arbitrary r and M can also 
be obtained from Eq.(25), and reads: 

oo 7 



e'^T^u^^^^ [ //r + l,//r + l,...,//r + l ] , [M_^ , A) = 



M times M times ( 53) 

00 V / 



In spite of their apparent complexity the l.h.s of the above equations can be 
straightforwardly handled by computer algebra systems [32] . 



6. Conclusions and outlook 



We have found exact analytical expressions for the generalized Stirling numbers 
and generalized Bell polynomials which appear in the normal ordering of powers of 
Laguerre-type derivative operators, and have provided a complete set of hypergeometric 
generating functions for these quantities. The combinatorial aspect of the problem 
was demonstrated by flnding an exact mapping between the normal ordering and an 
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enumeration of increasingly labelled, multivertex forests constructed according to a two- 
parameter (r, M) prescription. In this way analytical, numerical and combinatorial 
facets of this problem have been given a very complete treatment. We have also 
used generalized Dobihski relations to investigate the properties of these Laguerre-type 
differential operators. We provided a large number of operational formulas involving 
functions of Laguerre derivatives, which can alternatively be applied using the boson 
language. The framework developed above enables one to construct and analyze new 
coherent states relevant to nonlinear quantum optics, which will be the subject of 
forthcoming research. 
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8. Appendix: Sheffer-type operators 



We derive Eq.(35) with the help of methods developed in Ref.[8]. First, observe that 
L)(r, 1) = a'^a^^^ + ra^ from which it follows that D\r^ 1) = {a'^Y^^a + r{a'^y is an 
operator of Sheffer-type: D\r^ 1) = v{a'^) + q{a'^)a with q{x) = x^^^ and v{x) = rx^ . 
The normally ordered form of exp(Ai^'''(r, 1)) is obtained by solving the linear differential 
equations (Eqs.(2) and (3) of Ref.[8]) for T(A,x) and g{\^x) yielding 



and 



(1- Ar(at)-)iA 



1 



1 - Ar(at)^ * 

According to Eq.(29) of [8] the normally ordered form of e^^^^'-^^ is 



,AD(r,l) _ 



,ADt(r,l) 



= : i/(A,a)e^^ 



,at(T(A,a)-a) . 



which gives Eq.(35). 



[Al] 
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Figure 1. Building blocks (in the inset) and the associated graphs of order n = 1, 2 for 
Laguerre derivatives: (a) H = 1^(1, 1) = a^a^ + a and (b) H = D{2, 1) = a^a^ + 2a^. 




